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[57] ABSTRACT 

A method and apparatus for testing the operation of a 
complex stabilization circuit in a closed-loop system 
(12, 38) is comprised of a programmed analog (60-72) 
or digital (80-86) computing system for implementing 
the transfer function of a load (12), thereby providing a 
predictable load. The digital computing system employs 
a table stored in a microprocessor (84) in which pre- 
computed values of the load transfer function are stored 
for values of input signal from the stabilization circuit 
(38) over the range of interest. This technique may be 
used not only for isolating faults in the stabilization 
circuit (38), but also for analyzing a fault in a faulty load 
by so varying parameters of the computing system as to 
simulate operation of the actual load with the fault. 

7 Claims, 6 Drawing Figures 
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METHOD AND APPARATUS FOR TRANSFER 
FUNCTION SIMULATOR FOR TESTING 
COMPLEX SYSTEMS 

5 

ORIGIN OF INVENTION 

The invention described herein was made in the per- 
formance of work under a NASA contract and is sub- 
ject to the provisions of Section 305 of the National 
Aeronautics and Space Act of 1958, Public Law 85-568 10 
(72 Stat. 435; 42 USC 2457). 

FIELD OF THE INVENTION 

This invention relates to a method and apparatus for 
testing complex closed-loop systems by simulating the 15 
transfer function of a load in the system. 

BACKGROUND OF THE INVENTION 

While designing complex electronic systems, it is 
convenient and sometimes necessary to provide an 20 
adaptive load for testing. An example of such a complex 
electronic system is described by Michael J. Kavaya, 
Robert T. Menzies and Uri P. Oppenheim in a paper 
titled “Optogalvanic Stabilization and Offset Turning of 
a Carbon Dioxide Waveguide Laser,” published in the 25 
IEEE Journal of Quantum Electronics, Vol. QE-18, 

No. 1, January, 1982 at pages 19-21. That paper de- 
scribes a feedback loop employing the change in the 
discharge impedance of a CO 2 laser, as its output power 
is varied, to stabilize the frequency of the laser. A zero 30 
offset feature in the feedback loop allows continuous 
tuning of the stabilized frequency over a 300-400 MHz 
range. 

There are, of course, many closed-loop stabilization 
techniques for stabilizing the output frequency of a 35 
laser, such as a spectrophone stabilized laser, with line 
center offset, disclosed in U.S. Pat. No. 4,434,490, filed 
Mar. 31, 1982. It should therefore be understood that 
the method and apparatus for simulating a transfer func- 
tion of a load in testing complex systems disclosed 40 
herein is not limited to laser frequency control systems, 
although the system of that application will be used as 
an example. 

SUMMARY OF THE INVENTION 45 

In accordance with the present invention, the transfer 
function of a load to be stabilized by a feedback signal in 
a closed loop with a stabilization circuit is simulated 
while testing the operation of the stabilization circuit in 
the closed-loop system to verify its proper operation, 50 
and the limits of its proper operation may be determined 
by varying the parameters of the simulated load. 

The simulated load may be devised as an analog or 
digital computing means. In the latter case, software 
establishes the transfer function of the simulated load, 55 
which may therefore be more flexible for more complex 
loads. An equation is first written for the voltage output 
of the actual load as a function of an input voltage (out- 
put voltage of the feedback loop under test) and any 
reference voltages that may be applied. Digital comput- 60 
ing means are then programmed, such as with a look-up . 
table, to provide an output voltage, Vout, for every 
value of input voltage, V m from the stabilization cir- 
cuit. The input voltage from the stabilization circuit is 
converted to digital form and entered into the digital 65 
computing means for generation of an output voltage 
value as the transfer function of the load being simu- 
lated. The generated output voltage value is converted 
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from digital to analog form and applied to the stabiliza- 
tion circuit. This digital form of implementation is uni- 
versally applicable to testing any closed loop system for 
which the transfer function may be predetermined over 
a range of interest. 

The analog form of implementation is also based 
upon the known transfer function of the particular ac- 
tual load expressed by an equation. Analog computing 
means are then programmed to compute the output 
voltage of the test load using any reference voltages 
necessary in the computation. This output voltage is 
then applied to the stabilization circuit in the feedback 
loop of the complex system under test. 

The novel features that are considered characteristic 
of this invention are set forth with particularity in the 
appended claims. The invention will best be understood 
from the following example described with reference to 
the accompanying drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 is a schematic drawing of a spectrophone 
stabilized laser with line center offset frequency con- 
trol. 

FIG. 2 illustrates operation of the spectrophone stabi- 
lized laser of FIG. 1. 

FIG. 3 shows a plot of the first derivative of the 
Lorentzian line response of a spectrophone for various 
values of modulating amplitude. 

FIG. 4 shows a plot of the second derivative of the 
Lorentzian line response of a spectrophone for various 
values of modulation amplitude. 

FIG. 5 illustrates an analog computing means for 
generating the transfer function of the actual load for 
the stabilization feedback loop of FIG. 1 to test the 
operation of that feedback loop. 

FIG. 6 illustrates a digital computing means for gen- 
erating the transfer function of the actual load for the 
stabilization feedback loop of FIG. 1. 

DESCRIPTION OF PREFERRED 
EMBODIMENTS 

In order to present a complete description of the 
invention in a specific embodiment, the spectrophone 
stabilized laser with line center offset frequency control 
disclosed in the aforementioned application will first be 
described. 

Referring to FIG. 1 of the drawings, a CO 2 continu- 
ous wave (cw) laser 10 is shown as an example of an 
actual load 12 to be simulated in accordance with the 
present invention. A portion of the laser beam is di- 
rected through a spectrophone 14 by a beam splitter 16. 
The spectrophone is typically filled at a low pressure 
with the same gas as the laser. The gas pressure in the 
spectrophone determines the width of its absorption line 
which is a factor in the tuning control; the lower the 
pressure, the narrower the line. The width of the line 
may be defined through the parameter jl, which is the 
half-width at half maximum (HWHM), the frequency 
spread from peak center to a point at one half the ampli- 
tude of the peak. 

If the light beam through the spectrophone were 
modulated on and off as by a chopper, the electrical 
signal from a microphone 18 of the spectrophone 14 
would be proportional to the height of the original 
absorption curve. However, in accordance with the 
invention, a cw laser is modulated by a sinusoidal dither 
voltage applied to a piezoelectric transducer 20, such as 
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a lead zirconate-lead titanate (PZT) ceramic transducer, 
from an oscillator 22 applied through an adder 24 which 
sums two other voltages, one from a potentiometer 26 
providing a dc bias voltage Bl, and the other a correc- 
tion signal from a lock-in amplifier 28 superimposed on 5 
a second dc bias voltage B2 from a potentiometer 30 
through an adder 32 and integrated by an electronic 
circuit 34. The integrator is coupled to the adder by an 
amplifier 36. 

The voltage Bl is an offset voltage which is adjusted 10 
to tune the laser to the desired operating region within 
the spectral linewidth of the gas. The oscillator pro- 
vides the sinusoidal dither signal which is amplified by 
a circuit 40, the gain of which is adjustable, and applied 
to the PZT. This sinusoidal voltage dithers the length of 15 
the laser cavity ever so slightly, which in turn dithers 
the output of the laser correspondingly in a periodic 
fashion. In effect, the frequency is being dithered as 
shown along the absorption curve of graph A in FIG. 2. 

If the dithered spectrophone signal is fed into the 20 
lock-in amplifier and detected at a frequency gj, which 
is the original dither frequency, then the correction 
signal output Ol will be proportional to the first deriva- 
tive of the absorption curve, as shown in graph B of 
FIG. 2. The output signal, Ol, in this case will be zero 25 
at the center, positive to the left, and negative to the 
right. If desired, graphs B and C can be inverted with 
respect to the horizontal axis with amplifiers. If detec- 
tion is done at frequency 2 a>, the lock-in detection oc- 
curs at twice the dither frequency; the output signal Oi 30 
is proportional to the negative of the second derivative 
of the curve of graph A, or the negative of the deriva- 
tive of the first derivative in graph B, as shown in graph 
C. 

The first derivative has a zero-crossing at line center 35 
(point Pi) and two maxima at points P 2 and P 3 . The 
second derivative curve has a maximum at line center 
point P 4 and two zero crossings at points P 5 and Pe; the 
zero crossings are a certain distance from center corre- 
sponding to the maxima of the first derivative which has 40 
only one zero crossing that is always on line center. 

In previous stabilization circuits, the output of the 
lock-in amplifier was integrated directly, and the inte- 
grated signal was added into the PZT voltage, creating 
a compensation voltage at a stable operating point. The 45 
stable operating point is where the error signal E (FIG. 

1 ) is zero. Therefore, in detecting the dither signal using 
the first derivative, the only stable operating point is at 
Pi (graph B) where the output of the lock-in amplifier, 

O Ly is zero. If the second derivative technique is used, 50 
the only two stable operating points are at P 5 and P 6 
(graph C) where the output of the lock-in amplifier is 
zero. 

In the stabilization circuit 38, the zero-crossing points 
Pi, P 5 and Pa are not the only frequencies at which the 55 
system can be stabilized. If it is desired to stabilize at any 
given frequency (point on the detection curve of graph 
A), a bias voltage B2 (which is either positive or nega- 
tive) is generated using potentiometer 30 and added to 
the output of the lock-in amplifier 28 before integrating 60 
it as an error signal, E. 

The condition for stability is still E=0, but now for E 
to be zero, it is required that the lock-in amplifier output 
be equal to — B2 so that the sum (0£,-{-B2) be zero If it 
is thus required that the lock-in amplifier output be 65 
equal to — B2, and if B2 is set properly, it can be seen 
that on the first derivative curve of graph B, the stabi- 
lized frequency can be ranged anywhere between the 
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positive and negative peaks P 2 and P 3 with the adjust- 
ment of B2. If B2 is set at a level Q shown in FIG. 2, 
graph B, then the lock-in output has to be at level R, the 
negative of B2, which means that the laser operating 
frequency will be at the operating point T (graph A) 
and the laser will be stabilized there. 

In the second derivative technique, the initial operat- 
ing point may be, for example, at point P 6 . The operat- 
ing point will move either up or down the curve of 
graph C (FIG. 2 ) by adjusting voltage B 2 , since again 
the stability condition is that Ol~ — B2. 

The following is a summary of the operating regions 
in the first derivative technique. 


INTERVALS OF STABLE FREQUENCY OFFSET 
coDetection (1st Derivative) 

Lock-In Phase — 0 = 0* (Assume laser frequency increases 
as PZT voltage increases; and A3> 
A4, A5 > 0) 

POINT PI FOR B2 = 0 (as in the prior art) 

INTERVAL (1) FOR B2 < 0 
INTERVAL (2) FOR B2 > 0 

If the gain A3 of amplifier 36 (FIG. 1) is < 0, add 180° to 0. 


In this first derivative (a) detection) method, the lock-in 
amplifier would be set to 0 degrees (in phase) because 
the signal it detects at co would be in phase with the 
oscillator output A Sin coT, and stabilization in interval 
(1) if voltage B2 is less than 0, because that means the 
lock-in voltage is greater than zero, and stabilization in 
interval (2) if voltage B2 is greater than zero, because 
lock-in must be the opposite of it, as noted above and 
shown in FIG. 2 under the graph B. 

The following is a summary of the operating regions 
in the second derivative technique: 


INTERVALS OF STABLE FREQUENCY OFFSET 
2o) Detection (2nd Derivative) 

Lock-In Phase = 0 = 90° (same assumptions as in first 

derivative case above) 

INTERVAL (4) 

Subinterval (4A) if B2 > 0 
Subinterval (4B) if B2 < 0 
POINT P6 FOR B2 = 0 (as in the prior art) 

Lock-In Phase — 0 = 270° 

INTERVAL (3) 

Subinterval (3A) if B2 > 0 
Subinterval (3B) if B2 < 0 
POINT P5 FOR B2 = 0 (as in the prior art) 


In the second derivative technique there are four 
subintervals 3A, 3B, 4A, 4B of stabilization. The offset 
voltage B2 can be + or — , which means that the oper- 
ating point must be below or above the horizontal axis, 
respectively, but also in interval (3), the lock-in ampli- 
fier signal is rising in amplitude with increasing fre- 
quency and, in interval (4), the signal is falling in ampli- 
tude with increasing frequency. 

Not only does it matter in the second derivative case 
whether you are above or below the horizontal axis, 
which determines whether the offset voltage is — or -}-, 
but it also matters whether the signal is increasing or 
decreasing with increasing frequency; therefore, an- 
other parameter to consider is the phase of the lock-in 
amplifier. The 2 a> component of the spectrophone sig- 
nal into the lock-in amplifier is 90° out of phase with the 
oscillator output; therefore the lock-in amplifier is at a 
90° setting initially to stabilize in subintervals 4A and 
4B. Over this region the slope is the same as in the first 
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derivative case, failing with increasing frequency. If the 
lock-in amplifier phase is changed to 270°, stabilization 
is possible in subintervals 3A and 3B where the slope is 
opposite that of the first derivative case. 

Another important feature of the stabilization circuit 5 
38 is that the amplitude of the dither will determine how 
broad the detected curves are. As the dither amplitude 
is increased, the curves of graphs B and C in FIG. 2 get 
broader. What that means is that if it is desired to lock 
a little bit further out from line center on the frequency 10 
axis, the dithering amplitude can be increased to 
broaden the detected signal curves of graphs B and C. It 
is then possible to lock (stabilize) at that frequency. 

The spectrophone signal into the lock-in amplifier is 
proportional to the power of the laser. A technique for 15 
making that signal independent of the laser power uti- 
lizes another beam splitter 42 directing a portion of the 
laser beam into a power detector 44, obtaining the volt- 
age signal P l which is proportional to laser power, and 
then dividing the spectrophone signal S by that voltage 20 
signal in an analog dividing circuit 46, so that the output 
signal, V out, from the actual load is independent of the 
laser power. Thus, because the spectrophone signal S is 
proportional to laser power, after dividing by the laser 
power, the dependence is removed so that the depen- 25 
dence of the absorption curve, rather than the depen- 
dence of the output power of the laser with offset, is 
available in the output signal V out— &\S/ A 2 P l into 
the lock-in amplifier 28, where Ai and A 2 are fixed (or 
adjusted) gain factors of amplifiers 48 and 5D. 30 

It should be noted that when the laser is locked and 
stabilized at some point off line center, the amount of 
offset from center can be adjusted by three separate 
parameters: the voltage B2 offset; the amplitude of the 
dither, because that will widen or narrow the detection 35 
curve; and the pressure of the gas in the spectrophone, 
because that too will widen or narrow the absorption 
curve, thereby adjusting the amount of offset corre- 
spondingly. In practice, it will not be convenient to 
adjust pressure, once the instrument is set up for an 40 
experiment. Instead the pressure will be selected for 
some nominal breadth of the absorption curve. The 
breadth of the detection curve can then be modified 
during operation by adjusting the amplitude of the 
dither added by adjustment of the gain of the amplifier 45 
40. Note that coupling amplifiers 40, 48, 50 and 52, like 
amplifier 36, each have adjustable gain which may be 
unity, or may be greater or less than unity, and that 
amplifiers 48 and 50 may have their relative gains ad- 
justed to effectively set the range of the signal input to 50 
the lock-in amplifier. 

It will also be noted that the signal obtained with the 
first derivative technique is smaller than would be ob- 
tained with a light chopper; the signal obtained with the 
second derivative technique is smaller yet. This assumes 55 
that m < 1 where we define the peak-to-peak amplitude 
of the dither signal in units of laser frequency as 2m yx. 
These facts emphasize the importance of the use of a 
spectrophone as a sensitive transducer which makes 
these techniques possible. Also, the ability to choose the 60 
width of the absorption line through the choice of the . 
gas pressure in the spectrophone, allows greater flexibil- 
ity in the overall stability or “tightness” of the feedback 
loop. The technique is applicable to all three possible 
absorption lineshapes of the spectrophone gas fill: Lo- 65 
rentzian, Doppler, or Voigt. The modulation of an ab- 
sorption lineshape for the first two possibilities are de- 
scribed in the aforesaid patent application. 
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THEORY OF DERIVATIVE DETECTION 

Consider a general experimental configuration where 
the response or signal is affected by some experimental 
parameter and the shape of this effect is the Lorentzian 
lineshape as in the CO 2 cw laser of FIG. 1. Using fre- 
quency, v, as the experimental parameter, although it 
may be any parameter (e.g. voltage, magnetic field, 
etc.), the specific equations of the Lorentzian lineshape, 
its first three derivatives and their maxima are given in 
Table 1 shown in Appendix A. 

Assuming a small sinusoidal dither, at a modulation 
frequency of co, of the experimental parameter v, the 
resultant frequency is given by: 

v— vo+vi +myL sin o>t, (1) 

where vo is the Lorentzian linecenter frequency, v\ is 
the offset frequency from linecenter, yx, is the half- 
width at half-maximum (HWHM), and m is the modula- 
tion coefficient. The peak-to-peak modulation is 2myx . 
Inserting equation (1) into the Lorentzian equation 
yields the following equation: 


cr(v\,m,t) — - 


7 L 

(vj + myxsinci)?) 2 + 7l 

SyL 


( 2 ) 


it(v} -f yj) 


1 


1 + 


2vimy£sintD? + (myi) 2 s\n 2 o)t 

v 1 + ri 


for m<l, equation (2) is expanded using: 

= 1 - x + x 2 - x 3 + * 4 - . . . (x 2 < 1), 

to obtain 


cr(v\,m,t) = 


Sy L 


■ 


ri) 


1 + my i 


(4) 


— 2visintuf 
V 1 + 7L 


(myt) 2 


-1 

2<v? + yj) 


2*1 

(yj + y 2 D 2 


cos2 oit 
2(v? + yj) 


2V|COs2g)J 

O'? + yjj 1 


+ (myi ) 3 


3v]sino)f 

(v? + ri ) 2 


6vjsinoi/ 

(„1 + yfy 


visin3o>/ 2visin3ct)/ 

( V 1 + y]) 2 + (v 2 + yi) 3 

Equation (4) reveals that the signal consists of a dc 
component, a component at the modulation frequency 
co, and components at all the harmonics of co. The odd 
harmonics of co, including co itself, are in-phase with the 
initial dither signal, while the even harmonics are 90° 
out of phase. For values of m much less than 1, the high 
powers of m may be neglected. The signals at co and 2co 
become: 


+ 0 (m 4 ) 
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7T 


V\ 

[v\ 2 + 7L 2 ] 2 


myisinot + CKm?) 


and 


(5) 
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From the foregoing, it can be concluded that the 
Lorentzian lineshape for the actual load of FIG. 1 is 
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<T L {V) = 


_5 

7T 


7L 

(v - v 0 ) 2 + r/. 2 


(ID 


7 r 


(3vi 2 - yz. 2 ) 

!>1 2 + y £ 2]3 


( — m 2 y £ 2 c os2ct)/) 
4 


+ 0(m 4 ) 


( 6 ) 
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If the sinusoidal amplitudes of equations (5) and (6) are 
compared with the entries of Table 1 for v— vq— v\, it is 
evident that: 


in-phase signal at co 
first derivative 


= my L , 


and 


quadrature signal at 2a) 
second derivative 


—m 2 yi} 

4 - 


(7) 


15 


( 8 ) 


20 


Thus the o) signal is proportional to the first derivative 
and to myL and the 2o) signal is proportional to the 25 
second derivative and to (myL) 2 - Each signal may be 
increased by increasing m. However, this will decrease 
the resolution in v\ and eventually violate the assump- 
tion m<< 1. 

Equations (5) through (8) are approximations based 30 
on small dither amplitudes compared to the HWHM of 
the Lorentzian curve. However, the exact expression in 
equation (2) can be used with a Fourier series analysis to 
find the Fourier coefficients of the &> and 2co signal 
components. First, the offset frequency is expressed in 35 
terms of the HWHM by setting v\ The desired 

coefficients are then found to be: 


b\ 


25 

TiryL 



sin(<o/) 

1 + 03 + /wsiiuoO 2 


du 


a 2 = 


25 

TiryL 


T_ 

} cgsgff') du 

d-T I + 03 + /nsino>/) 2 
2 


(9) 

40 

( 10 ) 

45 


The integrals can be evaluated numerically. FIGS. 3 
and 4 show plots of yLbi/S and yLa 2 /S vs. £ for various 
values of the modulation coefficient m. The specific 50 
values of m are 0.1, 0.5, 1.0, and 2.0, where m = 2.0 
corresponds to the largest signal, etc. In addition to an 
increasing signal with larger values of m, it is evident 
that the curves widen in the horizontal or frequency 
offset dimension. (An analysis for large values of m 55 
shows that the o> signal is maximum when m = 2.2 for 
both the Lorentzian and Doppler line shapes.) From 
Table 1, it is seen that the peak values of the first deriva- 
tive (zero-crossings of the second derivative) occur at 
( v — vo) — ±yL/ v3 = ±0.577yL. These frequency offset 60 
coordinates are marked on the figures. Thus first or 
second derivative data should be used to estimate y l 
only for very small values of m, the modulation coeffici- 
ent. 

If large values of m are required because of signal-to- 65 
noise considerations, it may still be possible to derive 
yL by measuring the derivative signals versus frequency 
offset for two or more values of m. 


The essence of the present invention is in the method 
and means for providing an adaptive load for testing 
complex electronic closed-loop systems, such as a laser 
and frequency control circuit wherein the actual load is 
replaced by a transfer function simulator in which the 
transfer functions of the actual load are generated, using 
either analog or digital techniques. This simulation ena- 
bles feedback control loops under development to work 
into an adaptive load of known parameter values, in- 
stead of into an actual load that may itself have faults, 
making the loop difficult to analyze. 

Although the essence of the invention will now be 
described with reference to the CO2 cw laser of FIG. 1 
providing a spectrophone feedback signal as an example 
of an actual load for the particular stabilization circuit 
shown, it is recognized that there are many closed-loop 
stabilization circuits, either in use or proposed, for stabi- 
lizing the output frequency of a laser, and that the pres- 
ent invention has utility in the development of other 
complex feedback control systems. Consequently, it is 
intended that the claims be interpreted to cover this and 
other closed-loop stabilization systems in general. 

Referring now to FIG. 5, an analog simulator for a 
CO2 cw laser is disclosed for use as a test load in the 
development of the stabilization circuit 38, or to simu- 
late any actual load that exhibits a Lorentzian transfer 
function, for example. 

For the Lorentzian line transfer function ctl(v) given 
by equation (11), the output voltage is 


Vqut * 


TL 

Win - Yref) 2 + yi} 


( 12 ) 


The three inputs to this analog circuit shown in FIG. 5 
are: V/yy, the output voltage of the closed-loop stabiliza- 
tion circuit; Vref> the center voltage of the Lorentzian 
transfer function, an adjustable parameter; and yL, an 
adjustable voltage parameter which defines the width 
of the spectrophone absorption line. The output volt- 
age, V out, of the analog circuit is directed to the input 
of the stabilization circuit 38, which has the value AjS- 
/A2PL in the actual load, as shown in FIG. 1. 

The analog circuit shown in FIG. 5 is comprised of 
an attenuator 60, which brings the input voltage level 
into the range of an adder 62 where a reference voltage, 
—V ref, is added to the attenuated input, V'/M to form 
a difference signal V'in— Vref- This added reference 
voltage is adjusted to center the lineshape (the curve of 
graph A in FIG. 2). The output of the adder 62 is then 
squared by a multiplier 64 and added with the output of 
a multiplier 68 which squares the parameter yL that is 
adjusted for the desired transfer function width. The 
squared signal, yL 2 , is added with the difference signal 
(V'iN— V ref) in an adder 70. An output divider 72 then 
forms the output signal, V out • The analog signal thus 
functions in acordance with equation (12) for a voltage 
input, V/iv, from a stabilization circuit under test of the 
form 


Vjn=B[+A Sin ot -M 3 / E(t)dt 


03) 
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Note in the example of FIG. 1 that indeed only two 
voltages, V/wand V out, connect the stabilization cir- 
cuit to the actual load to be simulated. Consequently, it 
can be readily appreciated that the present invention 
can be applied to simulate any load in a closed loop 5 
having a stabilization circuit that receives a feedback 
signal and generates an output control signal for the 
load. As noted hereinbefore, the advantage of thus sim- 
ulating a load is that the stabilization circuit can be 
tested in the closed loop circuit without concern about 10 
the load itself having faults. 

Other laser configurations possible with the stabiliza- 
tion circuit 38 of FIG. 1 that can be simulated by imple- 
mentation of transfer functions include tuning the laser 
frequency by monitoring power transmitted through a 15 
gas cell. In either configuration (using a spectrophone 
or a gas cell) Stark electrodes can be provided on the 
monitoring device, and, instead of applying a modula- 
tion (dither) signal to the PZT voltage being corrected, 
it can be applied to the Stark voltage. Still another laser 20 
configuration for tuning frequency consists of monitor- 
ing reflection of the laser beam from a second optical 
resonator. Then either the laser PZT voltage or the 
resonator PZT voltage can be corrected and dithered, 
or else the laser PZT can be dithered while the resona- 25 
tor PZT is corrected, or vice versa. Yet another tuned 
laser configuration has already been referred to, namely 
tuning the laser frequency by monitoring voltage drop 
across the laser discharge (optogalvanic effect). In each 
case, the transfer function of the actual load is analyzed 30 
and then simulated, either in an analog circuit as for the 
example shown in FIG. 5 for the actual load 12 in FIG. 

1, or in a digital circuit shown in FIG. 6. 

Referring now to FIG. 6, the input voltage, V/jv, to 
the digital load simulator is first attenuated, as before, in 35 
a suitable analog attenuator 80 to bring the voltage 
within the range of operation of an analog-to-digital 
converter 82. Assuming a conversion to an eight bit 
binary number, the analog input signal may be digitized 
with a sensitivity of one part in 2 8 . Using equation (12) 40 


10 

cycle may require too much time for the desired fre- 
quency response of the closed loop. Consequently, the 
preferred implementation of the digital computation 
means would be to precompute each of the 2 8 values of 
Eot/rand store the values in a look-up table. For each 
value of V/n (attenuated as shown), there is thus a 
unique value immediately available for the output. This 
digital output is then converted to analog form by a 
digital-to-analog converter 86. This digital approach to 
simulation of the load has an advantage over the analog 
approach in that the load transfer function may be 
readily changed by simply changing the stored table. 
Thus a digital implementation provides a more flexible 
embodiment of the invention. It is universally applica- 
ble to testing any system for which the transfer function 
look-up table can be computed and entered into the 
table. 

In summary, there has been disclosed, with reference N 
to a specific example, a method and apparatus for elec- 
tronic simulation of a particular load in a closed loop 
characterized by predictable mathematical transfer 
functions. It is these functions, each of which may have 
its own set of parameters (Loren tzian, Doppler, Voigt, 
Airy, and the like), that the electronic computing means 
simulates, thereby providing a predictable load for test- 
ing a complex loop stabilization circuit which interacts 
with the load in the closed loop system. This technique 
is particularly useful in testing a circuit for stabilizing a 
load in a closed loop, such as a tuned laser, but may also 
be used for testing an actual load having a fault. For that 
purpose, one or more parameters of the load simulator 
are adjusted to match the fault of the actual load. In that 
way the fault of the actual load may be analyzed to 
determine what is required to correct it. 

Although particular embodiments of the invention 
have been described and illustrated herein, it is recog- 
nized that variations and equivalents may readily occur 
to those skilled in the art. Consequently, it is intended 
that the claims be interpreted to cover such variations 
and equivalents. 
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to compute V out, the voltage output of the simulated 
load for predetermined parameters V ref and jl, there 
is a unique value determined for each of 2 8 levels of 65 
VlN. 

The computation could be carried out on line using a 
programmed microprocessor 84, but the computation 


What is claimed is: 

1. A method for testing the operation of a stabiliza- 
tion circuit in a closed-loop system, where the stabiliza- 
tion circuit is designed to interact with an actual load in 
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response to an output signal, Vout, by generating a 
control signal, Win, from said output signal in accor- 
dance with the transfer function of said actual load, 
comprising the steps of determining the transfer func- 
tion of said actual load, programming computing means 5 
to implement said transfer function, utilizing said pro- 
grammed computing means in place of said actual load 
in said closed loop, and testing the operation of said 
stabilization circuit in said closed loop, thus determining 10 
the operation of the stabilization circuit in a closed loop 
system free of any faults that may be present in said 
actual load to be stabilized. 

2. A method as defined in claim 1 for analyzing a fault 

in an actual load wherein said transfer function imple- 15 
mented by said computing means includes at least one 
parameter, said method including the step of adjusting 
said parameter in said computing means to alter the 
characteristics of said load in order to emulate a faulty ^ 
operation of said closed-loop system with an actual load 
having a fault, thereby to determine the nature of said 
fault. 

3. Apparatus for testing the operation of a stabiliza- 
tion circuit having frequency stabilizing parameters in a 25 
closed-loop system for stabilizing a load, where said 
stabilization circuit interacts with an actual load in re- 
sponse to an actual load output signal, V out, to provide 

a feedback signal, Win, comprising a load simulating 
signal generating means adapted to be substituted for 30 
said load and supplied with said feedback signal and at 
least one load parameter, and operative to transmit an 
output signal, Wqut, as a function of said feedback 
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signal and load parameter for simulating said actual 
load, and means for adjusting said load parameter. 

4. Apparatus as defined in claim 3 wherein said load 
simulating signal generating means is comprised of an 
analog computing means responsive to said feedback 
signal and load parameter. 

5. Apparatus as defined in claim 3 wherein said load 
simulating signal generating means is comprised of an 
analog-to-digital converter for periodically sampling 
said feedback signal and converting it to digital form, 
digital computing means responsive to said feedback 
signal in digital form for computing said output signal in 
digital form as a function of said parameter, and digital- 
to-analog converting means for converting to analog 
form each output signal value computed in digital form. 

6. Apparatus as defined in claim 5 wherein said digital 
computing means is comprised of a look-up table in 
which precomputed values of said output signal are 
stored. 

7. Apparatus as defined in claim 3 wherein said actual 
load is comprised of a spectro phone stabilized laser 
having a Lorentzian absorption line response, and said 
computing means computes said output signal in accor- 
dance with the equation 


(V/N - v R ef ) 1 + yi? 

where yL is a parameter which defines the width of the 
spectrophone absorption line and V ref is the center 
voltage of the Lorentzian transfer function of said spec- 
trophone stabilized laser. 
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